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ABSTRACT 

\ We present the results of the asteroseismic analysis of the red-giant star KIC 4351319 

<^> ■ (TYC 3124-914-1), observed for 30 days in short-cadence mode with the Kepler satel- 

lite. The analysis has allowed us to determine the large and small frequency separa- 
tions, Ai'o = 24.6 ± 0.2 //Hz and 8v$2 = 2.2 ± 0.3 //Hz respectively, and the frequency 
of maximum oscillation power, y max = 386.5 ± 4.0 /iHz. 

The high signal-to-noise ratio of the observations allowed us to identify 25 independent 
p\ . pulsation modes whose frequencies range approximately from 300 to 500 /iHz. 

5-j ' The observed oscillation frequencies together with the accurate determination of 

the atmospheric parameters (effective temperature, gravity and metallicity) , provided 
by additional ground-based spectroscopic observations, enabled us to theoretically 
interpret the observed oscillation spectrum. 

KIC 4351319 appears to oscillate with a well defined solar-type p-modes pattern due 
to radial acoustic modes and non-radial nearly pure p modes. In addition, several 
non-radial mixed modes have been identified. 

Theoretical models well reproduce the observed oscillation frequencies and indicate 
that this star, located at the base of the ascending red-giant branch, is in the hydrogen- 
shell burning phase, with a mass of ~ 1.3 M©, a radius of ~ 3.4 R Q and an age of ~ 
5.6 Gyr. The main parameters of this star have been determined with an unprecedent 
level of precision for a red-giant star, with uncertainties of 2 % for mass, 7% for age, 
1 % for radius, and 4 % for luminosity. 

Key words: stars: individual: KIC 4351319 - stars: variables: solar-type - stars: 
oscillations - stars: red giants 



1 INTRODUCTION 

Oscillations excited by turbulent convection - solar-like os- 
cillations - have been successfully identified in cool main- 
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sequence and post-main-sequence stars with convective en- 
velopes. In the red giants, solar-like oscillati ons have been 
first detected by ground-based observa tions ()Merlinelll999l : 
iFrandsen et~al] |2002|; Jstello et alj l2004h and b y the WIRE 
satellite l|BuzasJl2000l ; lRetter et alj|2003l , 120041 ) . 

Difficulties in the identification of the ob- 
serve d modes, led sev eral aut hors jGuenther et all 

200d: iDziembowski et all l200ll : iTeixeira et all 120031 : 



Christensen- Dalsgaardl l2004al ~ to speculate on the pulsa- 
tional properties of the red-giant stars. Red giants are 
characterized by a deep convective envelope and a small 
degenerate helium core. Since the density in the core of 
these stars is quite large, the buoyancy frequency can reach 
very large values in the central part. In these conditions, 
g modes of high frequencies can propagate and might 
eventually interact with p modes giving rise to modes of 
mixed character. As a consequence, the spectrum of the red 
giants is predicted to have a quite complicated appearance, 
with a sequence of peaks uniformly spaced in frequency due 
to acoustic modes, and other peaks with less clear pattern 
due to mixed modes. 

The detection of solar-like oscillations in red giants 
has been well estab l ished by the space m ission MOST 
(|Barban et al.l 120071 ; iKallinger et~aH l2008al lbh. while the 
C0R0T satellite was able to de tect for the first time non- 
radial modes in red-giant stars (|De Ridder et alll2009T ) and 
to find solar-like oscillations in a very large sample of G 
and K giant stars ( Hekker et al 1l2009l : IKallinger" et al . 20101); 



iMosser et al. 1 120101 ) mainly lying in the core-helium-burning 
evolutionary phase (|Miglio et alj|2009l ). 

The high- q uality observation s of the Kepler mission 
(jBoruckil l2010l : iKoch et all I2010T ) enabled to extend the 
detection of solar-like oscillations from the red clump 
to the lower luminosity region of the red-giant branch 
jBedding et~aT]|2010l: iHekker et all l20ld : IHuber et all l20ld ; 
IKallinger et alfcoiOal ). where stars are still burning H in the 
shell. This finding enables us to study the structure and the 
evolution of the stars from the main-sequence phase up to 
advanced evolutionar y stages by means of asteroseismology 
jGilliland et al.ll2010l ). 

In this article we present the results of the analysis 
of the oscillation spectrum of the red giant KIC 4351319 
(TYC 3124-914-1) - nicknamed 'Pooh' (jMimell 1926b within 
the authorship team - observed by the Kepler satellite for 
30 days in short-cad ence mode (integration time ~1 min, 
iGilliland et~al1l2010D . Here, we also consider the theoretical 
interpretation of the observations and derive the asteroseis- 
mic estimates of age, mass and radius and of other structural 
characteristics, making use of the accurate atmospheric pa- 
rameters obtained by ground-based spectroscopy performed 
at McDonald Observatory. 

The star selected for this study is of particular in- 
terest because its oscillation spectrum shows an excess of 
power centered at freque ncy ~ 380 /jHz, higher than is typi- 
cal for Kepler red giants dBedding et al.l 20101; IHekker et ail 
|2010| ; IHuber et alll20ld : IKallinger et alj|2010al ). This indi- 
cates that the star is well below the red clum p in luminos- 
ity an d is still ascending the red-giant branch (|Miglio et al.l 
120091 ). Its relatively high frequencies require Kepler's short- 
cadence mode (see Section 4), which means that only rela- 
tively few stars with these characteristics are being observed 
by Kepler. KIC 4351319 therefore provides a very important 



opportunity to study this evolutionary phase and hence it 
has been selected for long-term observations by Kepler. 



2 ATMOSPHERIC PARAMETERS FROM 
SPECTROSCOPY 

Spectroscopic observations of KIC 4351319 in the visible 
spectral range have been carried out on 2010 July 26 
(HJD = 2 455 404.7756) with the Cross-Dispersed Echelle 
Spectrograph at the 2.7-m "Harlan J. Smith" telescope at 
the McDonald Observatory, Texas, USA, in the framework 
of the ground -based observational supp ort of the Kepler 
space mission (|Uvtterhoeven et al.ll2010a! lbT). 
A stellar spectrum, wavelength-calibrated and normalized to 
the continuum, was obtained using standard data-reduction 
procedures for spectroscopic observations within the NOAO 
/IRAF package. The resulting signal-to-noise ratio was 
SNR ~ 100. The spectral resolution as measured from the 
Th-Ar emission lines is about R = 65 000. 

In order to derive the effective temperature and surface 
gravity for our target, we minimized the difference between 
observed and synthetic profiles. As a goodness-of-fit test we 
used the parameter: 



2 

X = 



SIo 



where N is the total number of independent points in the 
spectrum, 7 b s and J t h are the intensities of the observed 
and computed profiles, respectively, and SI b s is the photon 
noise. The synthetic spectra were generated in three steps: 
first we compu ted the stellar atmosphere model by using the 
ATLAS9 code (jKuruczll 19931 ) . then the stellar spectrum was 
synthesized using SYNTHE (jKurucz fc Avrett I fl98lT ) and 
finally the instrumental and rotational convolutions were 
applied. The ATLAS9 code includes the metal opacity by 
means of distribution functions (ODF) that are tabulated 
for multiples of the solar metallicity and for various micro- 
turbulence velocities. 

We selected a set of vanadium and iron lines in the range 
between 6180 A and 6280 A, which are free of blending and 
whose atomic parameters are well kno wn and already ex- 
plored for temperature calibrations (see lBiazzo et al~ll2007l . 
and reference therein). For surface gravity we used the non- 
saturated line Cal at 6162.17 3 A, whose wings are very sen- 
sitive to gravity changes (see lGrevlll992f ). 

First of all, we computed the wsini of t he star by us- 
ing t he parameters from the KIC catalogue (jLatham et al.l 
l2005h reported in Table [T] We used SYNTHE to reproduce 
the observed metal lines, and our best match was achieved 
convolving the computed lines with a stellar rotational pro- 
file having usini = 6±l km s _1 , while the microturbulent 
velocity £ = lkms -1 was ob tained by using the r elatio n 
£ = £(T cff ,logg) published by I Allende Prieto et all (|2004T ). 
Then, with an iterative procedure starting with the values 
of the KIC cat alogue of T e ff and log fl and the solar abun- 
dances given bv lAsplund et all < |2009h . we obtained the best 
fit values reported in Table [1] Uncertainties in T e g and logg 
have been estimated within a 1 a level of confidence, as the 
variation in the param eters which increases the \ 2 by one 
(|Lampton et al. Ill976l ). The uncertainty in the iron abun- 
dance is the standard deviation of the weighted average of 



Solar-like oscillations from the depths of the red-giant star KIC 4351319 3 



Table 1. Basic parameters of KIC 4351319 





KIC and Simbad catalogues 


Present analysis 


V 


10.19 a 




v sin i 




6± lkms- 1 


T e ft 


4725 ± 200 K b 


4700 ± 50 K 


log g 


3.353 ± 0.5 dex 6 


3.3 ±0.1 dex 


[Fe/H] 


-0.563 ±0.5 dex b 


0.23 ±0.15 dex 



a http:/ /sirnbad. u-strasbg.fr/simbad/ 

b |http://archive. st sci.edu/kepler/kepler_fov/scarch.php 



Table 2. Summary of the inferred abundances in the atmosphere 
of KIC 4351319 



El 


[X/H] 


Na 


0.39±0.15 


Mg 


-0.02±0.14 


Si 


0.58±0.09 


Ca 


0.10±0.10 


Sc 


0.22±0.14 


Ti 


0.09±0.19 


V 


0.13±0.14 


Cr 


0.39±0.14 


Mn 


-0.01±0.14 


Fc 


0.23±0.15 


Co 


0.31±0.12 


Ni 


0.26±0.25 


Zn 


0.24±0.16 


Y 


-0.18±0.16 


Zr 


-0.18±0.15 


Ba 


-0.06±0.17 


La 


0.06±0.15 



the abundances derived from each single line. After having 
fixed Teff, logp and [Fe/H], we used SYNTHE to derive the 
abundance pattern of our target by the method of spectral 
synthesis. The results are summarized in Table 03 and plot- 
ted in Fig.[T] We then conclude that, according to estimated 
errors, the metallicity is slightly over the solar ones. In par- 
ticular, iron abundance is equal to [Fe/H] = 0.23 ± 0.15 dex. 



3 SOLAR-LIKE PROPERTIES IN 
RED-GIANTS STARS 

3.1 Large and small separations 

It is well known that the properties of solar-like oscillations 
in main-sequence stars can be describ ed by adopting the 
asymptotic development (Tassoul 1980), which predicts that 
the oscillations frequencies z/ nj ; of acoustic modes, charac- 
terized by radial order n, and harmonic degree I, for / <C n 
should satisfy the following approximation: 



Av 



I 

U+ 2 



(1) 



where e is a function of frequency and depends on the prop- 
erties of the surface layers and Av, known as the large fre- 
quency separation, is the inverse of the sound travel time 




Atomic Number 



Figure 1. Abundance pattern computed for all metals detected 
in the spectrum recorded with the "Harlan J. Smith" telescope at 
the McDonald Observatory. The horizontal dashed lin e represents 
the solar abundance as given in lAsplund et al ] 42009h . 



across the stellar diameter and is proportional to the square 
root of the mean density. 

Thus, the solar-like oscillations spectra should show a 
series of equally spaced peaks separated by Av between p 
modes of same degree and adjacent n: 



Av ~ AV\ = l/n+l, l - Vn.l- 



(2) 



In addition, the solar-like power spectra are character- 
ized by another series of peaks, whose separation is 5vi t i + 2, 
known as the small separation: 



8vi,l+2 = - v. 



(3) 



which is sensitive to the chemical composition gradient in 
the central regions of the star and hence to its evolutionary 
state. 

The determination of the large and small frequency sep- 
arations can provide asteroseismic inferences on the mass 
and the age of main-sequence and post-main-sequence solar- 
type st ars, using the so-called seism ic diagnostic C-D di- 
agram |Christensen-Dalsgaardl Il98&t ). However, these two 
seismic indicators alone do not represent an unambiguous 
diagnostic tool for more evolved stars, for which the rela- 
tion between the sma l l and large separation is almost linear 
(iBedding et all l20ld : iMontalban et al.1 l20ld : iMosser et all 
l2010h . On the other hand, it has been demonstrated that 
the stellar mass and radius of a solar-like star can suc- 
cessfully be derived, w ithin 7 % and 3 % respectively (e.g. 
iKallinger et al.ll2010bh . from measurements of the large sep- 
aration and Umax, the frequency at which the oscillation sig- 
nal reaches a maximum, b y usi ng the scaling laws g i ven b y 
iKieldsen fc Bedding] (|l995l ) and lBedding fc Kieldsenl (|2003h . 

Hence, in order to characterize in details the structure 
of a red-giant star, it is necessary to consider individual fre- 
quencies of oscillation and other oscillation properties, as we 
will discuss in the following. 
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3.2 The potential of the mixed modes 

The properties of solar-like oscillations are expected to 
change as the stellar structure evolves. According to Eq. fT} 
and considering that Au oc R~ 3 ^ 2 , oscillation frequencies of 
a given harmonic degree should decrease as the star evolves 
and the radius increases and should be almost uniformly 
spaced by Av at each stage of evolution. However, in sub- 
giants and red giants the radial modes seem to follow Eq. fl]) 
closely, but the frequencies of some non-radial modes appear 
to be shifted from the regular spacin g due to the occurrence 
of the so-called 'avoided crossing' (|Christensen-Dalsgaardl 
2004a). As the star evolves away from the main sequence, the 
core contracts and the radius expands, causing an increase 
of the local gravitational acceleration and of the gradients 
in the hydrogen abundance, and hence of the buoyancy fre- 
quency in the deep interior of the star. As a consequence g 
modes with high frequencies are allowed to propagate and 
can interact with p modes of similar frequency and same 
harmonic degree, giving rise to modes with mixed charac- 
ter, which behave as g modes in the interior and p modes 
in the outer envelope fAizenman et al.l 1 19771 ). The inter- 
action can be explained as the coupling of two oscillators 
of similar frequencies. The effect of the coupling becomes 
much weaker for modes with higher harmonic degree, since 
in these cases the gravity waves are better trapped in the 
stellar interior and hence better se parated from the region of 
propagation of the acoustic waves (iDziembowski et al.|[200ll ; 
IChristensen-Dalsgaardll2004al: iDupret et al.|l2009l). 

It has been found by iMontalban et al.l J2OI0T ) and ob- 
servationally demonstrated by Huber et al. l|2010l ). that the 
scatter of Z = 1 modes caused by 'avoided crossing' decreases 
as the star goes up to the red-giant branch: as the luminos- 
ity increases and the core become denser, the Z = 1 acoustic 
modes are better trapped and the oscillation spectra become 
more regular. Once the star ignites He in the core, the core 
expands and the external convective zone becomes shallower 
which has the effect of increasing the probability of coupling 
between g and p m odes again. 

Very recently, iBeck et alj (|201ll ) have demonstrated 
that the quality of the Kepler observations gives the pos- 
sibility to measure the period spacings of mixed-modes with 
gravity-dominated character which, like pure gravity modes, 
penetrate deeply in the core allowing to study the density 
contrast between the core region and the convective enve- 
lope and, like p modes, have amplitu de at the surface high 
enough to be observed. In particular, iBedding et al.l (|201ll ) 
found that measurements of the the period spacings of the 
gravity-dominated mixed modes, permit to distinguish be- 
tween hydrogen- and helium-burning stages of evolution of 
the red giants. 

The occurrence of mixed modes is then a strong indi- 
cator of the evolutionary state of a red-giant star and the 
fitting of the observed modes with those calculated by the- 
oretical models can provide not only mass and radius but, 
with a good approximation, an estimate of the age of the 
star. 



3.3 Sharp features inside the star 

Sharp variations localized at a certain acoustic depth in 
the structure of pulsating stars produce a distinctive quasi- 



periodic signal in the frequencies of oscillations. The charac- 
teristics of such signal are related to the location and ther- 
modynamic properties of the layer where the sharp variation 
occurs. Sources of sharp variations are for example the bor- 
ders of convection zones and regions of rapid variation in 
the first adiabatic exponent Ti = (dlnp/d In jo) a d, where 
the derivative corresponds to an adiabatic change, such as 
the one that occurs in the region of the second ionization of 
helium. 

In the main-sequence stars the signals coming from different 
sharp feat ures in the interior, overlap generating a complex 
behaviou r dMazumdar fc AntiallioOlh . In red giants, a recent 
study bv lMiglio et al.l <|2010h demonstrated that the oscilla- 
tory signal of the frequencies observed is directly related 
only to the second helium ionization zone. 

Several attempts have been made in order to iso- 
late the generated oscillatory components from the fre- 
quencies of oscillations or from linear combinations 
of them such as large separations and second differ- 
ences. In principle, this approach can be applied to 
determine th e properties of the base of the convec- 
tive envelope (|Monteiro~ Christcnscn-Dalsg aard. Thompson! 



2004) and in par- 



l200d : iBallot. Turck-Chieze and Garcfal 
ticular to inf er the helium a b undance in the stel- 
lar e n velope |Lopes et al ] [1997|; [Monteiro fc Thorni^ 



1998 : | Perez Hernandez and Christensen-Dalsgaardl Il99 



Miglio et al.l 120031 : iBasu et all 120041 : iHoudek fc Goughl 
2007ft . 



4 KEPLER OBSERVATIONS AND ANALYSIS 
OF THE OSCILLATION SPECTRUM 

KIC 4351319 has been observed by the Kepler satellite 
ijKoch et al.ll2010h during the Q0, Ql and Q2 runs in long 
cadence mode (integration time of ~ 30 min) and during 
the Q3 run in short cadence mode (integration time of ~1 
min). The data that we describe here were obtained dur- 
ing 30 consecutive days, starting on 2009, September 18, 
which corresponds to the first month of observing quarter 3 
(Q3.1). The Kepler short-cadence observation mode corre- 
sponds to a Nyquist frequency of 8497 fjJiz. In the present 
analysis, we used the SOC corrected data obtained by the 
Kepler team after reducing and correcting the raw photo- 
metric data for slop es a nd discontinu i ties a s described in 
IJenkins et all |2010l ) and lGarcfa et al.l |201ll ). The effective 
duty cycle resulted to be ~ 97% and the spectral window 
does not show aliasing artifacts that could interfere with the 
identification of real frequencies. 

The top panel of Fig. [2] shows the cleaned light curve 
normalized to the mean value. 

The Fourier analysis of the data set has been performe d 
by using the PERIOD04 package (|Lenz fc Bregerj 120051 ). 
which allows to extract the individual frequencies from large 
multiperiodic time series - also containing gaps - and can 
perform multiple-frequency least-squares fit up to several 
hundred simultaneous sinusoidal components. The resulting 
power spectrum of KIC 4351319 shows a clear power excess 
in the frequency range 300 — 500 /iHz with regularly spaced 
peaks (bottom panel of Fig. [2]), while no relevant features 
appear between 550 /iHz and the Nyquist limit. We have 
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Figure 2. Top panehLight curve corrected for slopes and discontinuities and normalized to the mean value. Bottom panel: Fourier 
spectrum in the region of p mode excess bump; harmonic degree of the excited modes are also shown 



considered only peaks in the spectrum with a signal-to-noise 
ratio (SNR) in amplitude greater than 3. 

In order to confirm the results obtained, the analysis 
of the spectrum was also independently performed by five 
additional teams which adopted the following methods: 

a) the power spectrum was smoothed with a Gaussian 
with a FWHM of 0.5 pHz and the frequencies of the peaks 
with SNR > 3 (in amplitude) were measured. 

b) peak-bagging methods, substantially based on the fit 
of Lorentzian profiles to the power density spectrum, using: 

• Maximum Likelihood Estimators (MLE), (see Ap- 
pendix A for details) 

• Bayesian MCMC l)Gruberbauer et~al 1 l2009h . (see Ap- 
pendix B for details) 

• Envelope autocorrelati on function analysis 
(M osscr fc Appourchauxll2009l ) 

Table [3] reports the set of frequencies confirmed, within 
the errors, by at least two of the peak-bagging solutions. 



Frequencies labeled with (•) have been detected by all the 
teams. 

The values of the radial order for / = are n — 
12 — 18 as obtained b y adoptin g the scaling law p roposed by 
iMosser et all (|2011aT ) (see also lHuber et al.ll2010l ). 

Three possible excited modes with 2 ^ SNR ^ 3, 
namely 302.32 /^Hz, 387.12 /iHz and 407.41 ^Hz have been 
also found, but their validity must be carefully checked by 
additional future observations. 

In order to estimate the large frequency separation, we 
first computed the autocorrelation function of the power 
spectrum in the region of p-mode power excess. With such 
an initial guess value of the large separation, we built an 
echelle diagram and identified the I = ridge. Afterwards, 
we computed t he co mb-response function, as defined by 
iKieldsen etaP ij 19951 ). for each frequency reported in Ta- 
ble [3] and for the three additional frequencies with 2 ^ 
SNR ^ 3, searching for the maximum value of the comb 
response in the range (15 — 30) fiHz. The result, reported in 
Fig. |3j shows that many modes lie around the mean value of 
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Table 3. KIC 4351319: observed oscillation frequencies with er- 
rors, their SNR and harmonic degree. The symbol • denotes 
modes confirmed by all the teams. 



Freq (/tHz) 



SNR I 



329.16 


±0.24 




5.4 





353.62 


±0.20 


• 


5.4 





378.34 


±0.11 


• 


13.1 





402.43 


±0.07 


• 


5.5 





427.32 


±0.12 


• 


7.5 





476.86 


±0.54 




3.6 





317.38 


±0.19 


• 


1.1 


1 


323.00 


±0.20 




3.4 


1 


332.98 


±0.30 




3.2 


1 


340.67 


±0.09 


• 


5.6 


1 


345.26 


±0.66 




3.2 


1 


365.40 


±0.07 


• 


8.9 


1 


370.09 


±0.15 




3.5 


1 


390.50 


±0.09 


• 


6.1 


1 


397.89 


±0.20 




3.8 


1 


417.30 


±0.11 


• 


5.1 


1 


440.20 


±0.10 


• 


4.7 


1 


326.55 


±0.25 




4.0 


2 


351.32 


±0.15 


• 


7.2 


2 


376.03 


±0.19 


• 


5.0 


2 


400.56 


±0.08 


• 


9.4 


2 


423.71 


±0.40 




3.4 


2 


425.80 


±0.32 


• 


3.2 


2 


448.65 


±0.40 




3.4 


2 


449.98 


±0.50 




3.1 


2 



(Ai^o) = 24.6±0.2 /iHz computed for the radial modes, while 
the other frequencies shows a quite large scatter indicative 
of their mixed mode character. 

Fig. [4] shows the region from 250 to 650 /iHz of the 
power spectrum folded with a spacing equal to the mean 
value of the large frequency separation. Two sharp peaks, 
corresponding to modes with / = and I — 2, as well as 
a broader peak, corresponding to modes with 1 = 1, can 
be clearly seen. The final mode identification in terms of the 
degree I, based on both the alignment in the echelle diagram 
(see Fig. I13[) and the result of the peak bagging is reported 
in Table O A value of the mean small frequency separation 
(8va2) = 2.2 ± 0.3 /iHz has been also derived from each pairs 
of v n ,o and v n -l,2 reported in Table[3] 



4.1 Background and p-mode power excess 
modelling 

After converting the power spectrum (ppm 2 ) to power spec- 
tral density (ppm 2 /iHz -1 ), by multiplying the power by the 
effective length of the observing run - estimat ed as the recip- 
rocal of the area under the spectral window (|Bedding et alj 
l2005h - a model fitting of the stellar background due to stel- 
lar activity and granulation has been performed. In order 
to simultaneously model the stellar background as well as 
the p-mode power excess, we fitted the observed power den- 
sity spectrum by a su perposition of white noise, two semi- 
Lorentzian functions (|Harvev| [i985) and a Gaussian func- 
tion re presenting the power excess hump (sec Kallinger et al. 
l2010bl , for details): 
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Figure 3. Large frequency separation as determined from the 
comb- response analysis applied to the frequencies of Table|3] (filled 
symbols) and to three possible frequencies with 2 ^ SNR ^ 3 
(open symbols). The horizontal line shows the mean value of 
24.6 ± 0.2 /tHz computed for the 1 = modes of Table [3] only. 
The reported uncertainty is the standard deviation of the mean. 
Different symbols refer to / = (stars), I = 1 (squares) and 1 = 2 
(circles) modes, respectively, as derived by the alignment in the 
echelle diagram. Triangles refer to mixed modes (see text for more 
details). 
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Figure 4. The power spectrum of KIC 4351319 folded at the 
average large separation, showing clear peaks corresponding to 
I = 0, 1 and 2. The small frequency separation Suq2 is also shown. 



P{U) 



2 



-4, 



^ (1 + vBt 



+ De~ 



(4) 



where v is the frequency, Ai, Bi and c; are the amplitudes, 
the characteristic timescales and the slopes of the power 
laws; E is the white noise contribution; D, Vmax and a are 
the height, the central frequency and the width of the power 
excess hump, respectively. 

We adopted as initial guess for the white noise the 
average of the power spectrum in the frequency range 
1000 — 2000 /iHz, where the photon noise is expected to 
dominate over the other components. The guess values for 
the remaining parameters have been evaluate d by scaling 
from typical solar value s (|Aigrain et al 1 12004 iPalle et al.l 
1 19991 : iHuber et al.|[20 09) for active region and granulation 
timescales. We smoothed the raw power spectrum density 
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Figure 5. Power density spectrum of KIC 4351319. Light grey: power density spectrum; dark grey: the raw spectrum smoothed with 
a 4/iHz boxcar taking only the independent averages; thick solid line: global fit according to the model described in the text; dashed, 
dotted and dashed-dotted lines: the two semi-Loretzians and the white noise component of the global fit, respectively. 



by applying a 4 //Hz boxcar using the ind ependent averages 
only, as suggested bv lGarcfa et al.l (|2009l ). in order to per- 
form properly the fit by means of a weighted least-squares 
proced ure. The weights we adopted, as in iGarcfa et al.l 
l|2009h . have been determined by t he uncertainti e s on each 
independent 4 /xHz average. As in iMathur et all (|2010h we 
adopted the IDL least-squares MPF1T package (provided 
by Craig B. Markwardt, NASA/GSFcEf), implementation 
of the Levenberg-Marquardt fitting algorithm. The results 
are shown in Fig. [5] 

The fit returned, among other parameters, the max- 
imum amplitude of the smoothed excess power P ma x = 
94.99 ± 0.3 ppm 2 //iHz. The corresponding frequency is 
Vm&x = 386.5 ± 4.0 pHz, where the uncertainty is given by 
the resolution of the smoothed spectrum. This value and its 
uncertainty are in a very good agreement with what found 
by fitting the raw spectrum by the envelope auto correlation 
function a nalysis dMosser fc Appourchauxll2009f ). 
Following iKieldsen et al.l (|2008al )7 we evaluated the maxi- 
mum amplitude per radial mode by adopting the calibration 
factor c = 3.16, whose value depends on the spatial response 
of the observations to modes with different degree relative 
to radial modes, and obtained A max ~ (27.19 ± 0.16)ppm. 



http://cow.physics.wisc.edu/~craigm/idl/idl.html 



4.2 Global asteroseismic parameters 

In order to extract a rough estimate of the global param- 
eters of the star, to be adopted as guess values for the 
model computation, we ado pted the scaling laws prov i ded b y 
IKieldsen fc Beddind (|l995l ) and lBedding fc Kieldsenl (|2003l) , 
relating the observed mean large frequency separation 24. 6± 
0.2 fiHz and i/ m „ = 386.5 ± 4.0 /xHz to the fundamental pa- 
rameters of the star. We obtained: M/M© ~ 1.35 ±0.09 and 
R/Kq ~ 3.44±0.08. We also derived the expected radial or- 
der of the mode with maximum amplitude in the spectrum, 
n max ~ 14. This latter value is in agreement with the ra- 
dial order identification of the radial mode with higher SNR 
reported in Sect. 4. 



5 ON THE CHARACTERIZATION OF THE 
STRUCTURE OF KIC 4351319 

5.1 Evolutionary models 

Given the identified pulsation frequencies and the basic pho- 
tospheric parameters, we faced the theoretical challenge to 
interpret the observed oscillation modes by constructing 
stellar models which satisfy the observational constraints. 

We assumed the effective temperature and gravity as 
calculated in Section 2 respectively, T c ff = 4700 ± 50K, 
logg = 3.3 ± 0.1 dex (TableJlJ. 
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Figure 6. Evolutionary tracks plotted in the H-R diagrams (the 
upper panel shows T c ff — logg, the lower panel shows T c ff — L) 
calculated for increasing values of the mass, while all the other 
parameters are fixed. The mctallicity is Z = 0.03, the initial hy- 
drogen abundance Xq = 0.7, and a = 1.8. Here, no additional ef- 
fects are included. The rectangle defines the two-sigma error box 
for the observed gravity and effective temperature. In the bottom 
panel, black dashed lines indicate the uncertainty in the observed 
effective temperature while observed values of luminosity are not 
known. 



We calculated several grids of theoretical structure 
models for the star by using the ASTEC evolution code 
(IChristensen- Dalsgaard 20083) by varying the mass and the 
composition in order to match the atmospheric parameters 
available. 

All the models hav e been calculated with th e OPAL 
2005 equa tion of state dRogers fc Navvonovl I20021 ). OPAL 



Rogers! 1 1996l\ , and the NACRE 



opacities 

clear reaction rates ( Angulo et al.l I1999T ) . Convection was 
treated according to the mixing-length formalism (MLT) 
|Bohm-Vitensd [l958|) and defined through the parameter 
a = I /H p , where H p is the pressure scale height and a is 
assumed to vary from 1.6 to 1.8. 

Inclusion of overshooting or diffusion outside the con- 
vective core during the main-sequence phase did not pro- 
duce any appreciable effect on the oscillation frequencies of 
models in such evolved phases of the evolution. 

A crucial input quantity is the iron abundance, whose 
logarithmic value relative to the solar one has been taken 
as [Fe/H]= (0.23 ± 0.15) dex, as determined in Sect. 2. 
The initial heavy element mass fraction Z can be calculated 



from the relation [Fe/H]= \og(Z/X) - \og(Z/X)®, where 
(Z/X) is the ratio at the stellar surface and t he solar value 
is (Z/X) e = 0.0245 (|Grevesse fc Noels!ll993h . Thus, it has 
been assumed Z/X — 0.04 ± 0.01. In the present modelling 
we neglect the difference in the relative abundances of the 
heavy elements between the star and the Sun (cf. Fig. 1); 
although this should be taken into account in future, more 
detailed investigations, the effects on the opacity and hence 
the model structure are likely to be contained within the 
assumed range of uncertainty in Z. The uncertainty in the 
observed value of Z introduces an uncertainty in the de- 
termination of the mass whose value, considering only the 
observed spectroscopic parameters, seems to be limited to 
the range M = (0.8 - 1.5) M Q . 

Extra mixing outside the convective region was ob- 
tained by adding turbu lent diffusion from t h e con vective en- 
velope, as described in lProfntt fc Michaudl (|l99ll ). by using 
the parameterized formulation: 



(p/Pb c J" 



(5) 



where pb cz is the density at the base of the convective en- 
velope and -Dmax is a parameter which sets the maximum 
diffusion. 

Additional evolutionary models were calculated by in- 
cluding overshoot beneath the convective envelope by a dis- 
tance £ ov = o v^p,b cz! where -ff p ,b cz is the pressure scale 
height at the base of the convective envelope and ct ov is a 
nondimensional parameter. 

The resulting evolutionary tracks are characterized by 
the mass M, the initial chemical composition and a mixing- 
length parameter. Fig.[6]shows a series of evolutionary tracks 
obtained for different masses and fixed initial composition 
and plotted in two H-R diagrams, representing respectively 
the effective temperature-gravity plane and the effective 
temperature-luminosity plane. 

The location of the star in the Hertzsprung-Russell di- 
agram identifies KIC 4351319 as being in the post-main- 
sequence phase of evolution at the beginning of the ascend- 
ing red-giant branch. It has a small, degenerate helium core, 
having exhausted its central hydrogen, and it is in the shell- 
hydrogen-burning phase, with a very deep convective zone, 
extending from the base located at about rb cz — 0.2i? to the 
photosphere. 

Figure[7]shows the internal content of hydrogen in three 
models of this star, which include respectively extra mixing 
by turbulent diffusion from the convective envelope, over- 
shooting from the convective envelope with a ov = 0.2 and 
without any extra mixing effect. The internal hydrogen pro- 
files show that the core, where the hydrogen is exhausted, is 
very small and a sharp variation, as a step-like function, 
occurs at the base of the convective envelope located at 
mb c , = 0.28M and rb cz = 0.23i? for the model with no 
extra effects considered. Models with inclusion of turbulent 
diffusion show a smooth profile of the abundance of the hy- 
drogen with no sharp variation at the base of the envelope. 

We expect to be able to distinguish among the different 
internal characteristics of the structure by studying the oscil- 
lation properties of this star. F urthermore, we wish to verify 
wheth er, as demonstrated by IChristensen-Dalsgaard et al.l 
(|l993h for the case of the Sun, the inclusion of diffusion in 
the models results in a better agreement between theory and 
observations. 
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Figure 8. Evolution of adiabatic frequencies with age and with effective temperature of a model computed with a mass M = 1.32Mq 
and Z = 0.03. The dotted lines correspond to modes of degree £ = 0, and the solid lines to modes with I = 1. The frequency and effective 
temperature ranges correspond to the observed ones. The grey dots indicate the acoustical cutoff frequency. The vertical line indicates 
Model 1 of Table H 



T eii (K) 

4800 4750 4700 




Age (Gyr) 

Figure 9. Evolution of adiabatic frequencies for modes with I = 2 as function of the age and the effective temperature for a model of 
KIC 4351319 computed with a mass M = 1.32M0 and Z = 0.03. Here solid lines correspond to modes with even radial order, grey dotted 
lines correspond to modes with odd radial order. The frequency and effective temperature ranges correspond to the observed ones. The 
vertical line indicates Model 1 of Table [4] 
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Figure 7. Hydrogen content in three models of KIC 4351319 
computed with M = 1.32M Q , X = 0.7, Z = 0.03, a = 1.8. 
The base of the convective envelope located at mt, cz = 0.28M 
and ru = 0.23R, is shown by the dot-dashed line. The solid 
line corresponds to Model 1 in Table|4]with no additional effects. 
Dotted line corresponds to Model 3 calculated with inclusion of 
strong turbulent diffusion (Z5 max = 40000cm 2 s^ 1 ). Dashed line 
corresponds to Model 4 of Table [4] calculated with overshooting 
from the convective envelope. 



Figure 10. Propagation diagram from the centre to r = 0.3 R 
for three models. The solid black line represents the buoyancy 
frequency N for the Model 1 in Table [4] without extra-mixing 
effects, the solid grey line represents the buoyancy frequency N for 
the Model 2 in Table [4] with the addition of a weak diffusion from 
the convective envelope, the dotted line represents the buoyancy 
frequency N for the Model 3 in Tablc[4]with an efficient turbulent 
diffusion from the convective envelope. The dashed lines represent 
the Lamb frequencies Si for 2 = 1,2, 3. The dot-dashed horizontal 
lines define the limits of the range of the observed frequencies. 



5.2 The seismic properties of the models 

In order to investigate the ob served solar-like oscillations, we 
used the ADIPLS package (jChristensen-Dalsgaardl l2008bl ) 
to compute adiabatic oscillation frequencies with degree 
I = 0,1,2, 3 for all the models satisfying the spectroscopic 
constraints. 

Figure [5] shows the evolution of frequencies computed 
for an evolutionary model of KIC 4351319 calculated with 
M — 1.32Mq, Z = 0.03, without additional extra mixing 
effects. The ranges in frequency and effective temperature 
have been chosen to correspond approximately to the ob- 
served ranges. The location of the acoustical cut-off fre- 
quency, decreasing with increasing age, at the top of the 
atmosphere in the model has been indicated by grey dots. 

As it has been already explained in Sect. 3, accord- 
ing to Eq. fT}, the plot should be characterized by frequen- 
cies which decrease as the star evolves and almost uniformly 
spaced by Av at each stage of evolut ion. In the case of more 
evol ved stars, as it has bee n shown bv lDi Mauro et all (120031 ) 
and iMetcalfe etal] (|2010h . while the radial modes seem to 
follow closely Eq. ([T}, the frequencies of the modes with 
I = 1 show the typical oscillating behaviour due to the pres- 
ence of the mixed modes. In addition, Fig. [8] shows that at 
each stage of the evolution the distance between two fre- 
quencies of adjacent orders become smaller and smaller as 
the frequency decreases: while the upper part of the panel 
hosts mixed modes with prevalent p-mode character, below 
350 ^Hz the panel appears more crowded by mixed modes 
with predominant g-mode character and the general trend 
of frequencies decreasing as the star evolves does no longer 
hold. 

In particular, pure g modes of high radial order are lo- 
cated at very low frequencies and high radial orders, below 
300 /iHz. 

Fig. [9] shows the evolution of the frequencies for I — 2 
modes as function of the age and the effective temperature. 



It shows that I = 2 modes, similarly to what happen for 
I — 1 modes, undergo avoided crossings during the evolu- 
tion and the theoretical oscillation spectrum is as a result 
very crowded. The presence of mixed modes with prevalent 
g mode character is also evident for 1 = 2: below 400 /iHz 
frequencies do not appear to decrease with the decrease of 
the effective temperature. 

The behaviour of the mixed modes inside the star can 
be understood by considering the propagation diagram in 
Fig. [10] obtained for three models of KIC 4351319, calcu- 
lated respectively without extra mixing effect and with the 
inclusion of a weak and an efficient turbulent diffusion be- 
low the convective envelope. The characteristics of the non- 
radial modes depend on the separation between gravity and 
acoustic domains, defined respectively by the Brunt- Vaisala, 
(TV) and Lamb frequencies (Si). The Lamb frequencies ap- 
pear close to the g-modes region. Thus, the gravity waves, 
trapped in a region not so well separated from the region 
of propagation of the acoustic waves, can interact with the 
p modes of similar frequencies. In addition, it is possible to 
notice that at nearly the same value of effective temperature 
and luminosity, the model without extra-mixing effect shows 
a second maximum in its buoyancy frequency located at the 
base of the convective envelope (r = 0.23i?). This maximum 
becomes smaller and totally disappears as the effect of the 
helium diffusion beneath the base of the convective envelope 
becomes stronger. The characteristics of this second maxi- 
mum in TV reflects the behaviour of the hydrogen abundance 
inside the star as seen in Fig. [7] and it will be discussed in 
Sect. 5.4. 

It is important to notice that, according to the theory, 
the spectra of g modes are denser than that of p modes: 
many g modes can interact with the same p mode of close 
frequency. 

In order to have an idea about the modes expected to be 
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Figure 11. Normalized inertia plotted against frequency for 
Model 1 described in Table [4] Black circles are employed for / = 
modes. Blue triangles connected by solid lines are employed for 
I = 1 modes. Green squares connected by dotted lines are em- 
ployed for I = 2 modes. Red diamonds connected by dashed lines 
indicate 1 = 3 modes. 



visible, amplitudes of all the calculated modes were roughly 
estimated by making the following considerations. The to- 
tal energy of a mode of radial order n and harmonic de- 
gree I can be expressed as £„i = A 2 „iE„i, where A„i is the 
surface amplitude, and E„i is the inertia of the mode. For 
stochastically excited modes it is generally expected that the 
mode energy is independent of the harmonic degree at fixed 
frequency. It follows that the amplitude of a mode can be 
estimated relative to a radial mode of the same frequency, 
such that the modes would have the same total energy (e.g. 
IChristensen-Dalsgaard et al.|[l99l iDupret et al.ll20091 ): 



where Aq and Eo have been obtained by interpolating to 
the frequencies of radial modes. Figure [TT] shows the nor- 
malized inertia for the frequencies of Model 1 of Table [4] of 
KIC 4351319. It shows that radial modes have a very low 
inertia relative to most of the nonradial modes and are the 
modes which will show highest amplitude at the surface. 
Some modes with I = 1, / = 2 and I = 3 have an inertia 
similar to that of the radial modes E n i ~ Eo, and appear 
quite regularly spaced in frequency. These modes which pre- 
serve strong p-mode character are solar-like oscillations and 
are predicted to be visible at the surface with a quite high 
amplitude. 

The effect of coupling, typically expected only for 1 = 1, 
strongly appear also for modes with I > 2. As the iner- 
tia grows, during the coupling the gravity character of the 
modes predominates on the acoustic character. Modes with 
the highest inertia are evanescent towards the surface and 
hence have the lowest probability to be observed. Thus, due 
to the low amplitude, the probability to detect any mixed 
modes for / = 3 appears quite low, but it is quite likely that 
we can observe many of the mixed modes for I = 1 and 
prob ably a few of the mixed modes for I = 2 dBeck et al.l 

In order to clearly distinguish the nature of the mixed 
modes, it is useful to study the behaviour of the displace- 



ment eigenfunctions of the modes. Figure[T2lshows the radial 
variation of the eigenfunctions of three mixed modes calcu- 
lated for Model 1 of KIC 4351319 described in TableU The 
upper panels of Fig. [12] show the eigenfunction for a mixed 
mode with high inertia and hence with a predominant grav- 
ity nature, characterized by a quite large amplitude in the 
core. The middle panels of Fig. [T2l show the eigenfunction for 
a mixed mode with both gravity and pressure character. The 
lower panels of Fig. fT2] show the eigenfunction for a nearly 
pure acoustic mode, with an inertia comparable to that of 
the radial modes: the amplitude in the core is negligible. 

5.3 Interpretation of the observed oscillation 
spectrum 

Among all the possible models, we selected a few models able 
to reproduce the set of all the observed frequencies of Table 
[3]and the observed values of the large and small separations. 
The characteristics of some models that satisfy the observed 
constraints are given in Table [4] According to the stellar 
evolution constraints, given the match with the observed 
oscillation properties, and with the use of all the possible 
values of mass and metallicity, our computations show that 
the age of KIC 4351319 is (5.6 ± 0.4) Gyr, with a mass M = 
1.30±0.03 Mq, a radius R = 3.37±0.03 H@ and a luminosity 
L = 5.1 ±0.2 L . 

The values of mass and radius which have been obtained 
by direct modelling of the individual frequencies are better 
constrained and well compatible with those obtained by the 
scaling laws as described in Sect. 4.2. 

A detailed comparison between the theoretical oscil- 
lation spectra for Model 1 and Model 5 of KIC 4351319 
reported in Table U and the observed data is provided in 
the echelle diagrams of Fig. 1131 The size of the symbols is 
proportional to the theoretical oscillation amplitudes of p 
modes, relative to the amplitudes of radial modes with the 
same frequency (see Eq. [BJ. The results show, as explained 
in previous sections, that the observed modes are I = pure 
acoustic modes, and / = 1 and / = 2 nearly pure p modes 
and some I = 1 and I = 2 g-p mixed modes. Some non-radial 
modes with mixed gravity-pressure character have an inertia 
so low as to propagate up to the surface and appear to be- 
have like solar-like oscillations. Very few mixed modes, with 
a quite high inertia, keep their gravity character, although 
the combination with a p mode enhance their amplitude so 
that they can be observed at the surface and in the echelle 
diagram appear to depart strongly from the regular solar- 
like pattern. 

Pure g modes at lower frequencies, have not been detected 
with the present observations. Solar-like oscillations with 
/ = 3, although theoretically predicted, have not been de- 
tected, probably due to geometrical reasons which produce 
the cancellation of the observed signal. 

We found that there is a very good agreement between 
observed and theoretical frequencies. In particular we are 
able to reproduce all the I = 1 mixed modes. The presence 
of mixed-modes with / = 2 needs to be better investigated: 
unfortunately none of our best models, shown in Fig. 1131 
is able to reproduce the I = 2 mode with frequency v = 
423.71 ^iHz. In addition, the results show that it is quite 
difficult to determine which of Model 1 and Model 5 best 
reproduces the observations. Perhaps the ambiguity will be 
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Figure 12. Radial variations of cigcnfunctions for three selected mixed modes calculated for Model 1 in TableU The left panels show 
the behaviour of the eigenfunction from the centre to 0.1 R. The right panels show the eigenfunctions from 0.1 R to the surface. The 
harmonic degree I, the radial order n and the frequency v of the modes considered are indicated inside each panel. 

Table 4. Relevant parameters, the mass M, the Age, the luminosity L, the effective temperature T e ff , the initial hydrogen abundance Xq, 
the surface metallicity Z s , the surface radius R, the location r c j-, of the base of the convective envelope in units of R, the mixing-length 
parameter a, the diffusion coefficient -D max and the overshooting parameter ctov for a set of models of KIC 4351319. 
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solved by the identification of modes at lower frequencies 
and in particular if the reality of the mode of 302.32 /iHz 
(see Sect. 4) will be confirmed or not. 

Another important issue is related to the possibility of 
observing rotational splittings in this star. We estimate that 
the observed rotational velocity v sin i ~ 6 km s~ 1 , quite high 
for this kind of objects, might produce a rotational splitting 
of ~ 0.8 /iHz for the modes of I = 2, m = ±2 in case of rigid 
rotation. Such a rotational splitting is about twice the value 
of the frequency resolution of the 30 d data set. To check 
for the signature of rotation we have re-run the frequency 
determination now adding rotational split components for 
I = 1 and I — 2 modes. The rotation frequency and the in- 
clination angle were allowed to vary in the range between 



— 3 /iHz and 0° — 90° re spectively, following the approach 
of lGizon fc Solankil (120031 ). A comparison of the global like- 
lihoods of the two models (with and without rotation) in- 
dicates that the rotation-free model is about 70 times more 
likely than the model including rotation. 
Indeed, the presence of two modes closely spaced in fre- 
quency does not prove that we deal with rotational splitting 
although this possibility cannot be excluded. Moreover, the 
two pairs of close frequencies of I = 2 suspected for evi- 
dence of rotational splittings (448.65 fiRz with 449.98 ^Hz 
and 423.71 /iHz with 425.80 /iHz), appear to be not sepa- 
rated by a same amount in frequency, indicating that they 
cannot be both appearance of rotational splittings. 

We finally notice that in order to reproduce the ob- 
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Figure 13. Two echelle diagrams based on observed and computed frequencies, plotted with Ai/ = 24.6 /iHz. The upper panel shows 
results for Model 1 in Table [4] while the lower panel shows results for Model 5 in Table [4] The open symbols represent computed 
frequencies, while the filled symbols represent the observed frequencies of Table [3] Circles are used for modes with I = 0, triangles for 
I = 1, squares for I = 2, diamonds for I = 3. The size of the open symbols indicates the relative surface amplitude of oscillation of the 
modes. 



served modes, we did not need to c orrect the frequencies fo r 
the surface effects, as sugg ested bv lKieldsen et all (|2008bh . 
but the small departure of the I = modes at high fre- 
quency might be indicative of the need of such a correction. 
In fact, the structure of the near-surface regions of the stars 
is quite uncertain: there are still substantial ambiguity in 
modelling the convective flux, defining an appropriate equa- 
tion of state to describe the thermodynamic properties of the 
stellar structure, as well as in the treatment of non-adiabatic 



effects. The conditions and limits of the applicability of the 
surface effect are still unknown and will be the subject of 
future studies. 

Since there is no evidence for fast rotation, we exclude any 
possible effect on the frequencies due to magnetic activity. 

It was possible to determine the evolutionary state of 
this star through direct modelling, which places it in the 
hydrogen-shell-burning phase at the base of the ascending 
red-giant branch. This star is substantially less evolved than 
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Figure 15. First adiabatic exponent as function of the acoustic 
radius for Model 1 of Table [4] The region where helium under- 
goes its second ionization corresponds to the second minimum in 
Ti. The location of the base of the convective envelope is also 
indicated. 



the red giants, wh ose I = 1 grav i ty-mo de pe riod spacings 
were m easured by iBedding et al.l (|201lf ) and iMosser et al.l 
This is reflected in Fig. [8] where we see that the 
spacing of the g-dominated modes (the upward sloping fea- 
tures) is comparable to the p-mode spacing. This fact, to- 
gether with the very strong mode bumping, makes it is quite 
hard to see a regular spacing in the 1=1 modes of any sin- 
gle model. This makes it difficult to measure a clear g-modc 
period spacing. However, we see from the list of observed 
modes in Table 3 that the consecutive I = 1 modes have pe- 
riod spacings in the range 35 -50 s, which is c onsist ent with 
a hydrogen-shell-burning sta rlBedding et all (120111 ). Apply- 
ing the method described bv lMosser et al.l ~ 2011bl ) confirms 
a period spacing of about 40±4 s. 



5.4 Sharp features inside KIC 4351319 

Figure [14] shows that a weak oscillatory motion is present 
in the calculated frequencies for the acoustic modes which 
follow closely the asymptotic law, namely the radial modes 
and the non-radial modes with very low inertia. This oscilla- 
tory signal, as explained in Sect. 3.3, is produced only by the 
region of the second ionization of helium, which induces a 
local minimum in the first adiabatic coefficient Fi and hence 
a sharp variation in the sound speed. 

Figure 1151 which shows the behaviour of the first adi- 
abatic coefficient as function of the acoustic radius t = 
Jq 1 dr/c(r), indicates that in the models of KIC 4351319 the 
second helium ionization zone occurs at about tneii/T ~ 0.7, 
while the base of the convective envelope is located at 
tb C!lI /T ~ 0.08, too deeply to produce an effe ctive signal into 
the solar- like frequencies (|Miglio et alj2ufch . The properties 
of the helium ionization zone, once determined from the os- 
cillation frequencies, may be used to constrain the structure 
of the star, in particular the envelope helium abundance. 
Further studies and developments are required in order to 
determine how a different content of He and a different equa- 
tion of state for the computation of the models will result in 
a different oscillatory behaviour in the theoretical frequen- 



Secondly, we have also considered the possibility to iso- 
late the signal coming from the base of the convective en- 
velope, by looking at the observed modes which penetrate 
deeply inside the star, namely the m ixed modes with a high 
inertia. This can be studied following iMiglio et al.l (|2008t l. by 
comparing models modified by adding turbulent diffusion or 
overshooting at the base of the convective envelope. The re- 
sults have shown that the inclusion of turbulent diffusion, 
or overshooting, below the convective envelope produces a 
displacements only in the mixed modes for / = 1 and / = 2 
modes, while all the other frequencies are not modified as 
shown in Fig. 1161 Thus, although we have not observed high- 
order g modes, we can still use the mixed modes to probe 
qualitatively the interior of this star. However, it is impos- 
sible to distinguish the effect on the frequencies due to the 
inclusion of diffusion from that of overshooting. 



6 CONCLUSION 

KIC 4351319 has been observed for 30 days by the Kepler 
satellite. These observations have yielded a clear detection 
of 25 modes identified with harmonic degrees I — 0,1,2 
between 300 and 500 /iHz with a large separation A^o = 
24.6 ± 0.2 ^Hz and a small separation 8vo2 = 2.2 ± 0.3 pHz 
respectively. 

The oscillation spectrum of this star is characterized 
by the presence of a well defined solar-like oscillations pat- 
tern due to radial acoustic modes and non-radial nearly pure 
p modes equally spaced in frequencies. The nearly pure p 
modes have mixed gravity-pressure character with an iner- 
tia so low as to propagate up to the surface and appear to 
behave like solar-like oscillations, although they penetrate 
deeply to the core. 

The oscillation spectrum also showed evidence for the 
presence of mixed modes with strong gravity-mode charac- 
ter. Those modes sound the internal region of this star and, 
besides constraining the stellar age, carry tight information 
about the location of the convective envelope and the condi- 
tion in the core. This enabled to define characteristics of this 
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star with an accuracy which cannot be reached by the use 
of only the global asteroseismic parameters, Av and i/ max . 

Thus, in this paper we have addressed the problem of 
identifying the structural properties and the evolutionary 
state of KIC 4351319 by using the results of the analysis 
of the oscillation spectrum and the atmospheric parameters 
provided by supplementary ground-based spectroscopic ob- 
servations. 

Detailed modelling of this star, in the attempt to match 
both the asteroseismic and spectroscopic constraints, al- 
lowed us to determine the main parameters with an un- 
precedent level of precision for a red-giant star, with un- 
certainties of 2% for mass, 7% for age, 1 % for radius, and 
4 % for luminosity. In the end, we are able to conclude that 
this star is in the red-giant phase of the evolution, with a 
mass M = (1.30 ± 0.03) M©, an age of 5.6 ± 0.4 Gyr, a ra- 
dius R = 3.37 ± 0.03 R Q and a luminosity L = 5.1 ± 0.2 L Q . 
The uncertainties obtained for the stellar parameters come 
from the degeneracy of the models describing the stars in 
the red-giant phase of the evolution. 

Only the detection of g modes, which better sound the 
condition of the core, could improve the results obtained by 
the present analyses, while strong theoretical efforts need to 
be done in order to achieve a better description of the stellar 
structure. 

It is clear that KIC 4351319 represents an excellent can- 
didate for long-term observations, which have been already 
scheduled for Kepler future runs. This will allow us to open 
a window on the understanding of this phase of the evolu- 
tion and to unveil all the raised questions including those 
relative to the identification of the mixed modes, to the 
signature of the ionization of the He and to the presence 
of diffusion or overshooting below the convective envelope. 



In particular, the possibility to detect rotational splittings 
in KIC 4351319 appears very fascinating and represents an 
unique opportunity for sounding the internal dynamics of a 
red-giant star. 
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APPENDIX A: MAXIMUM LIKELIHOOD 
ESTIMATORS (MLE) 

According to this method, the p-mode parameters are esti- 
mated by finding the best fit between a modelled and the 
observational power d ensity spectrum by us ing a maximum- 
likelihood technique (|Anderson et alj Il990l ) . Two different 
models have been assumed: 

Al Lorentzian with convolution of the window 
function 

The fitting of the non-oversampled power spectrum has 
been done in two steps. First the background, for which 
we use a power law fit to account for long term effects, such 
as granulation, rotation or activity related to spots, is fit- 
ted. Then, we compu ted a global maximum likelihood fit 
(|Anderson et al.l ll990T l to all oscillation modes, where the 
final model M is a convolution of the model with the power 
spectrum of the observed window function of the data set, 
normalised to unit total area (see Eq. IA1|) . This takes the 
redistribution of power caused by gaps in the data into ac- 
count. 

M= V 77 — r- — — I * window (Al) 

\j * + ((k'cenj — V)/-Bj) 2 y ^ ' 

In this model the previously determined background is kept 
fixed, while each of the jth oscillation peak is fitted with 
the height (H), the central frequency (v ce n) and HWHM 
(B) and the noise level as free parameters. 

The selection of the oscillation modes is based on a sta- 
tistical test of the binned power spectrum. For this test we 
bin the power spectrum over intervals of three frequency 
bins. Then we compute the probability of the power in the 
binned power spectrum to be due to noise. This probability 
is computed using a \ 2 distribution in which we take the 
width and number of the bins into account in the degrees of 
freedom. Frequencies at which the probability of the power 
not being due to white noise is larger than 95% are selected 
as candidate oscillation frequencies. After performing the fit, 
we also verify that in the ratio of the observed to the fitted 
spectra no prominent peaks are left. Therefore we compute 
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the relative height for the investigated frequency range for 
which the probability of observing at least one spike with 
this height due to no ise is 10%, following the formulation by 
IChaplin et all i|2002l ) and references therein. 
This is also described in iHekker et all (|2010u . 

A2 Sum of Lorentzians or sine 

If the modes are resolved (i.e. the mode width is larger than 
the frequency resolution) we chose to model the oscillation 
spectrum by a sum of Lorentzians plus a background model 
that counts for the signal which is not due to the p modes 
such as the instrumental noise or stellar background. In the 
frequency range considered here, the background has been 
modelled by a straight line (i.e. a + bp). When the mode 
width is of the order of the frequency resolution we change 
the sum of Lorentzians by a sum of sine. 

The modelled power spectrum used to match the data 

is: 

Q 

P("k) =^2M(n,u k )+B(u k ), (A2) 

n=l 

where Q is the number of oscillation modes, n the radial 
order, Uk the Fourier frequencies, and B(v^) the background 
noise in the power spectrum, modelled as a straight line, i.e. 
a + buk with a and b, two constants and free parameters of 
the fit. -M(n, v k ) is a Lorentzian when the mode width is 
larger than the frequency resolution: 

M(n,u k ) = R n -i (A3) 

1 I l?(v-k-Vn)\ 



where H n is the height of the Lorentzian profile, v n is the 
oscillation mode frequency and F is the mode line width 
(FWHM) with T = 1/(ttt), r being the mode lifetime. 

When the modes are not resolved, A4(n,^k) is a sine 
function: 

M(n, u k ) = H n sinc 2 (nT(^ k - u n )) (A4) 

where T is the total length of the observing run. 

The power spectrum is fitted "globally" over a fre- 
quency range corresponding to the detected excess power. 
The free parameters of the fitting process are: 

• the height of the Lorentzian profile, H n , n being the 
radial order. A single height parameter is fitted per mode. 

• the frequency for each mode, v n . 

• the line width for each mode (FWHM) with F n = 
l/(nT n ), r n being the mode lifetime. A single width has been 
considered. 

• the parameters a, b describing the background model as 
mentioned above. 

The mode-parameter la error bars are derived from the 
Hessian matrix. No oversampling has been used in the com- 
putation of the power spectrum for the fitting procedure in 
order to minimize the correlation of the points. The degree 
of each mode is identified using echelle diagrams. 



we use the previously determined stellar activity and gran- 
ulation signal (Eq. 4 without the Gaussian) as a fixed back- 
ground. The visually identified mode profiles are parame- 
terized by their central frequencies, mode heights, and life- 
times. As it can be seen in Fig. 2, the mode profiles are quite 
narrow. To prevent the algorithm from over-fitting the data 
and to keep the number of free parameters to a reasonable 
amount we use two lifetime parameters, one for 1=0 and 2 
modes and one for 1=1 modes. Furthermore, we assume the 
individual mode heights to follow a Gaussian envelope pa- 
rameterized by a single central frequency and width but with 
individual heights for each mode degree. For the 25 modes in 
Tab. 3, this results in a total of 32 free parameters. For the 
fit we use a Bayesian MCMC algorithm (|Gruberbauer et alj 
l2009h that delivers probability density functions for all fitted 
parameters and their marginal distributions, from which we 
compute the most probable values and their la uncertain- 
ties. During the fit the mode frequency parameters are al- 
lowed to vary independently within ±2/iHz around the value 
inferred from a visual inspection of the spectrum. The life- 
times are kept between 1 and 50 days. Whereas the centre 
and width of the Gaussian envelope are kept within 0.9 and 
1.1 times i/ mM and a (Eq.4), respectively, the height for each 
mode degree is allowed to vary independently between and 
2 times the highest peak in power density spectrum. 



APPENDIX B: BAYESIAN MCMC METHOD 

The p-mode parameters are estimated by fitting a sequence 
of Lorentzian profiles to the power density spectrum, where 



